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Abstract 

In general the onset of hybrid inflation requires an extremely homogeneous field config- 
uration at the end of the Planck era. The field <j> orthogonal to the inflaton a must be zero 
with high accuracy over a range exceeding the initial horizon size ~ ni~j^ by about two orders 
of magnitude. We consider a supersymmetric model that permits the decay of the oscillating 
field 4> into light particles. We study the field evolution and find that the requirement of 
extreme homogeneity can be relaxed. However, the field <fi must still be smaller than the 
inflaton by a factor of order 1 over a region far exceeding the initial horizon size. 
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1 Initial conditions for hybrid inflation 



An attractive realization of the inflationary scenario is obtained when hybrid inflation |I] is 
embedded in the context of supersymmetry ||]— |J. Flat directions with non-zero potential 
energy density appear without fine-tuning and are not lifted by quantum corrections. The 
COBE observation of the cosmic microwave background anisotropy constrains the properties 
of the model along the inflationary trajectory [10]— |12j. On general grounds, one expects the 
inflationary energy scale V 1 ^ (determined by the vacuum energy density during inflation) to be 
at least two or three orders of magnitude smaller than the Planck scale F] [13]: 



V*/4/ e l/4 _ ? x 1Q 16 Q e y ( X ) 



Here e is a "slow-roll" parameter [12] that must be much smaller than 1 during inflation. 



The question of initial conditions for hybrid inflation has been addressed recently [ 14 ] — JTJ 
In general, the onset of inflation requires a region of space with a size of a few Hubble lengths 
where the fields take almost constant values, so that the gradient energy density is negligible 



compared to the potential energy density [17]. For hybrid inflation, the requirement is much 
stricter 

Let us consider a supersymmetric model of hybrid inflation |2[ || based on the renormalizable 
superpotential 

W = -fi 2 S + \S$ 2 . (2) 

It is invariant under a global i?-symmetry with the superfields S and $ carrying i?-charges 1 
and respectively. The scalar components of these superfields can be written in the form 

5 = 4=, $ = (3) 

Here we have used an .^-transformation in order to make S real. The potential is given by the 
expression 

V = M 4 - XfJ 2 (fx - 4) + ^ (fx + 4) 2 + A V (f x +<%). (4) 

The supersymmetric vacua are located at a = 0, <j)\ = 4fi 2 /\, <p2 = 0. The potential has 
a flat direction along the a axis, for a > <7j ns = /x/V% with the non-zero vacuum energy 
7(<r,0i = 4>2 = 0) = /i 4 that can support inflation. The flatness of the potential along the 
a axis is lifted by radiative corrections, resulting from the breaking of supersymmetry by the 
non-zero value of V(a, = 02 = 0). The precise value of these corrections is not important for 
our discussion. They induce the slow roll of the inflaton field once inflation has set in, but have 
a very small effect on the question of the onset of inflation. For this reason we neglect them in 
the following. 

We assume that a Robertson- Walker metric is a good approximation for the regions of space 
with uniform fields that we are considering. The evolution of the fields is given by the standard 
equations 

dV 

a + 3H&= - —, (5) 
dV 

01,2 + 3^l, 2 = -TT— , (6) 



x We refer to the "reduced" Planck scale m m = Mpi/V&k, M Pi = 1.22 x 10 19 GeV. 
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Near the onset of inflation, the fields 0x2 oscillate around zero with amplitudes [0i 2]- The 
equation of motion for a reads 

a + 3i^ = -A 2 ([0i] 2 + [0 2 ] 2 )a, (8) 

where we replaced 2 and 02 by their average values (0f 2 ) ~ \ \4>\,2^ during an oscillation. 

Two time scales characterize the solutions of this equation. The first one is related to the 
"friction" term and is given by = 3H/2. Near the onset of inflation, the Hubble parameter 
is dominated by the non-zero vacuum energy V = [i 4 , where the scale [i is constrained by the 
COBE observations to be two or three orders of magnitude below mp\. This gives 

t - 2 mpl (Q\ 

The other time scale is obtained if we neglect the "friction" term and consider the oscillations 
of the a field. One-fourth of the period is the typical time for the system to roll to the origin 
and away from an inflationary solution. It is given by 
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(10) 



Inflation sets in if t osc tn and the a field settles at a constant value on the a axis. This 
requirement can be written as 

yVi] 2 + [<fe] 2 ^ Vsir /x 

\i ~ 4 A m pi 

For a scale \x consistent with the COBE observations {jijmpi ~ 1CP 3 -1CP 2 ) and couplings of 
order 1, the value of the fields 0i 2 orthogonal to the inflaton must be zero with an accuracy of 
two or three orders of magnitude in units of fi. This condition of extreme homogeneity must be 
satisfied throughout a region of space with a size of the order of the Hubble length ~ mpif/j, 2 . 
Otherwise, the fields in different parts of this space region will evolve towards very different 
values. In one part they may settle in the valley along the inflaton axis, while in another they 
may end up at the minima of the potential. As a result, large inhomogeneities will appear at 



scales smaller than ~ H^ 1 that will prevent the onset of inflation [16]. 

On the other hand, the earliest time at which one could start considering regions of space with 
classical fields is when the Planck era ends and classical general relativity becomes applicable. 
The initial energy density is near m Pl . For a theory with couplings not much smaller that 1, 
the initial field values within each region, as well as the initial field fluctuations, are expected 
to be of order mp\. It is very difficult for such an initial field configuration to evolve into the 
extremely homogeneous configuration that is required at the onset of hybrid inflation. 

In fig. 1 we present the evolution of the fields <r, 0i for a theory described by the potential 
of eq. (U) with A = 0.3, fj,/rnpi = 10 -3 . The initial values of the fields have been chosen 
a/mpi = 3, <fii/mpi = 0.6, 02/w-p; = 0.5. They are near mpi, with 0i, 02 smaller than a, 
so that the evolution starts near the flat direction. The initial time derivatives of the fields 
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have been set to zero and the evolution equations ©-(0) have been integrated numerically. We 
observe an initial stage during which the a field is reduced, while 4>\ and oscillate rapidly 
around zero with quickly decaying amplitudes. This quick decay is due to the large initial value 
of H, induced by the large average value of V in eq. (|). When the fields become much smaller 
than their initial values, H drops significantly and the "friction" term in eqs. dBf), (P) becomes 



less important. The amplitudes of the <p± 2 fields are much larger than the bound of eq.(ll). As 
a result, a rolls below zero and into negative values. Eventually, the system of fields enters a 
stage of oscillations around zero with slowly decaying amplitudes. The frequency of oscillations 
of 4>i t 2 is larger than that of a because of the difference of their effective masses induced by the 
last term in the right-hand side of eq. (||). 

It is clear from the above discussion that the onset of hybrid inflation requires the fine-tuning 
of the initial field configuration in such a way that the amplitudes of ^12 are smaller than the 
bound of eq. ( |i~l| ) when the energy density has dropped from ~ m pi to ~ /x 4 . In general, this 
leads to a requirement of extreme homogeneity for the initial field configuration at the end of the 
Planck era (with an accuracy of up to ten orders of magnitude in units ofmpj)-an "unnatural" 
assumption Jl~6[ . 

What should be the initial scale of the homogeneity at the end of the Planck era? The size 
L of a homogeneous region evolves proportional to the scale factor and, in terms of the Hubble 
length, is given by 

J*. (12) 
L H \poJ 

The parameter w determines the relation between energy density and the mean value of the 
pressure (p = wp). For a system of massive oscillating fields, or a matter-dominated Universe, 
w = 0, and for a radiation-dominated Universe w = 1/3. At the onset of inflation, when 
p/rripi ~ p 4 /rripi ~ 10~ 8 -10~ 12 , we must have LH ~ 1. For this we must require LqHq ^ 10- 
100 at the end of the Planck era, when po/m pl ~ 1. This means that the initial extreme 
homogeneity must extend far beyond the typical size of regions that can be considered causally 
connected. 

A feature of the field evolution that was not taken into account in ref. Ilql concerns the 



possible decay of the oscillating fields 4>i,2 into light particles [18|. As these fields are coupled to 
the inflaton, they can produce a particles. However, the production rate through the coupling 
\ 2 a 2 (d>\ + (f) 2 ) is proportional to the square of the amplitude of the (f>\ t 2 oscillations and becomes 
negligible very quickly. As a result, this particle production does not alter the picture we sketched 
above. 

The possibility exists, however, that the oscillating fields 4>i t 2 produce other light particles 
with a large decay rate. In this letter we consider a supersymmetric model for which such a 
scenario is realized, and show that in this situation the requirement of the extreme homogeneity 
can be relaxed. Namely, provided that the preinflationary reheeting is effective, inflation can 
proceed from the space region in which after the Planck era the initial field values are <j)\ 2 ^ 



a ~ mpi. However, still the size of this region should be larger than the initial horizon size 



m Pi by about two orders of magnitude. 



2 Particle production before the onset of inflation 

We consider the renormalizable superpotential 

W = -p 2 S + \S<Z> 2 + nm 2 + m^ 2 . (13) 
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It is invariant under a global i?-symmetry with the superfields S, <I>, carrying i?-charges 1, 
0, 1/2 respectively. We have not included a term even though it is allowed by the R- 

symmetry, because it can be eliminated through a redefinition of the superfields. Moreover, 
this term is absent if the superfields $ and belong in non-trivial real representations of some 
symmetry group under which S is a singlet. We assume that the two mass scales [/,, m are 
comparable and the couplings A, k not much smaller than 1. 

The scalar components of the superfields can be written in the form 



(14) 



Here we have used an /^-transformation in order to make S real. The potential is given by the 
expression 



V = n - Xfi 2 



+ (^1 + ^2 



2m z 



X 2 

+ T 



+ 2\f2~Km<p\ 



+ 



+ nXa 



The supersymmetric vacua are located at a = 0, cf>\ ■ 
potential has a flat direction (independent of a) for <f>i 
direction the masses of the cj)\ 2 fields are given by 



0i (V? - 4>i 



V/a, h 

-- 02 = Ipl ■ 



+ 202^1^2 



(15) 



Mi 



Mi 



2X 2 a 2 



2X^, 
2X 2 a 2 + 2X^ 2 . 



0, ^1 = 0, ip 2 = 0. The 
■02 = 0. Along the flat 



(16) 



An instability appears for 



^ 2 / ^-2 



X 



;i7) 



which can trigger the growth of the (pi field. 

The coupling cjxj/ 2 permits the decay of the 0i 5 2 fields into ^ particles. If these fields 
perform coherent oscillations, they can decay into scalar or fermionic particles. In the most 
favourable situation, a is slowly varying in the region a S> Oi ns and the frequency of oscillations 
is \/2Xa ^> 2m. The decay rates of each one of the 0i t 2 fields into scalar or fermionic ^ particles 
are equal f\ and 



r = r s + r / , 



r 



1 



/ 



4\/27r 

The evolution of the fields is given by the equations 



k 2 Xu. 



a + 3H& 

H, 2 + (3if + r)0i 2 



dV 
fa' 
dV 

501,2 ' 



(18) 

(19) 
(20) 



In our model the fermionic ^ particles are Majorana spinors. 
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dv 2 r~5~ 



^1,2 + 3#Vl,2 = " ^ " « y ^2 ^1,2, (21) 

Pfl + 4ff^ = r((0f> + (^>), (22) 



tf 2 



(5)"=4(?^ + H- (23) 



The quantities (</>f 2 ) correspond to the average energy of oscillations of the <j)\2 fields over one 



cycle ]19| . The gas of particles with energy density pp consists of scalar and fermionic Vl/ particles 
with mass 2m. The second term in the right-hand side of eq. ( pi] ) results from the interaction 
of the ipip fields with the gas of scalar particles (see below). 



We have integrated numerically the above equations for the potential of eq. (15) with 
k = 0.5, A = 0.3, n/mpi = 10 -3 , m/mpi = 10~ 2 . The initial values of the fields have been 
chosen a/mpi = 3, 4>i/mpi = 0.6, 4>ijmpi = 0.5, i^i/mpi = 0.1, tp2/mpi = 0.2, They are near 
mpi, with (/>i, ip\, ip2 smaller than a, so that the evolution starts near the flat direction. The 
parameters of the a, 0i >2 system have been chosen equal to those of the previous section (see fig. 
1). The initial time derivatives of the fields have been taken equal to zero. The initial energy 
density is po/m pl ~ 0.512 and the initial value of the Hubble parameter Hq/tupi ~ 0.413. 

In fig. 2 we present the evolution of the field <pi. We observe that it settles quickly in a pattern 
of regular oscillations around zero. The frequency of oscillations is set by the mass in the first of 
eqs. (|l~6|). The amplitude of the oscillations decays quickly because of the large initial values of 
the Hubble parameter and the decay rate of eqs. (]I8|). After a time tmpi ~ 5 x 10 2 , the evolution 
of <p\ is influenced by its interaction with the oscillating ipi ^ fields through the last term in the 
third line of eq. (|l5|). As a result, <f>\ develops an oscillatory mode with characteristic period 
~ ir/m. In fig. 2 we also present the evolution of the a field. We observe that it rolls slowly 
towards the origin because of the effective mass generated by the first term in the third line of 
eq. (|i~5l). As the amplitude of the 4>i t 2 oscillations drops, the time derivative of a is reduced 
because of the "friction" term 3H& in its equation of motion. 

The oscillating (fri^ fields produce an equal number of scalar and fermionic particles. The 
scalar ^ particles have a large self-interaction rate and their gas is expected to reach thermal 
equilibrium quickly. However, the interactions of the fermionic particles involve additional 
heavy $ particles and their superpartners, and thermalization is not possible for them. The 
temperature T of the scalar ^ gas is given by 

PR 4 

with an effective number of degrees of freedom = 2. The interaction of the '01,2 fields with 
the thermal bath of the scalar \& particles is established through the term k 2 (i/j 2 + ^2) 2 /4 in 
the potential. It results in an effective thermal mass for these fields 



>»Un =y T2 = K V^2V^. (25) 



We have included this contribution in the equations of motion (^H), as it can be much larger 
than the zero-temperature mass 2m. The interactions of the tpi^ fields with the fermionic \P 
particles are negligle, as they involve heavy <]? particles and their superpartners which are not 
expected to be abundant. 
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In fig. 3 we present the evolution of all the fields during the time between tmpi = 10 3 and 
tmpi = 2 x 10 3 . After tmpi ~ 10 3 , the fields 0i2> ^12 follow a pattern of oscillations with 
characteristic frequency set by the mass of the ij)x,2 fields = M^ 2 = 2m = 0.02 mp\. The 
small oscillations of 4>i,2, with high characteristic frequency set by their mass in eqs. ([l6|), quickly 
die out. The amplitude of the low- frequency oscillations is reduced because of the "friction" term 
~ 3H in the equations of motion. The Hubble parameter is dominated by the energy density pp 
of the gas of ^ particles. At t\jmpi = 1.5 x 10 3 this energy density is [pp] 1 /m pl ~ 3.28 x 10~ 7 , 
while the (potential and kinetic) energy density of the fields is / m pi — 5.86 x 10~ 9 . The 
Hubble parameter is Hi/mpi ~ 3.34 x 10~ 4 . The field a has a value crx/mpi ~ 0.978 and a time 
derivative &i/m Pl ~ —7.26 x 10~ 5 . 

We have found numerically that, in the time interval between tmpi ~ 10 3 and tmpi ~ 
1.5 x 10 3 , the (kinetic and potential) energy density of the fields pf falls proportionally to a 
power 4 of the scale factor R. The reason is that the dynamics of the fields during this time is 
strongly affected by the quartic terms in the potential and there is still significant decay of the 
4>i : 2 fields into ^ particles. At times tmpi ^ 1.5 x 10 3 , pf falls proportionally to approximately 
the third power of R. This is due to the fact that the amplitude of the field oscillations becomes 
sufficiently small for the quartic terms to be irrelevant, while the particle production practically 
stops. The radiation energy density falls as R~ 4 until a time tmpi ~ 2 x 10 3 . At this point the 
typical energy of a scalar ^ particle E ~ T is expected to be comparable to its mass. Similar 
conlusions can be reached for the typical energy of a fermionic particle. Based on the above 
considerations, we assume that, after tmpi 1.5 x 10 3 , the Universe becomes matter dominated 
and the energy density of the particle gas falls as -R -3 , similarly to the field energy density. As a 
result, the ratio of the field energy density pf to the radiation energy density pp stays constant 
e / ~[p / ] 1 /[pi i ] 1 ~1.79xl0- 2 . 

The evolution of the a field after t\/mpi = 1.5 x 10 3 is very simple. The right-hand side of 
eq. ( |l~9|) is negligible compared to the "friction" term proportional to 3H in the left-hand side. 
The Hubble parameter as a function of time is given by H = 2/ [3(1 + w)t\. The parameter w 
determines the relation between energy density and the mean value of the pressure {p = wp). 
For a system of massless oscillating fields, or a radiation-dominated Universe, w = 1/3. For a 
system of massive oscillating fields, or a matter-dominated Universe (such as the case we are 
considering), w = 0. The solution of eq. (|i"9| ) is given by 

\+ W _g_ / -l+w \ 

a(t) = ax + (T!^ t{ + ™ [tr w ~ t—j ■ (26) 

Asymptotically, the a field is expected to settle down at a value aj/mpi = 0.869. We have 
verified numerically the evolution of a according to eq. (|2^) and that the Universe is always 
dominated by the energy density of the gas of ^ particles. 

Inflation is expected to start at a time tmpi ~ 10 6 , when pp ~ p 4 . The size of the amplitudes 
of the 4>i t 2 field oscillations at this point is crucial for the onset of inflation. They must satisfy 
eq. ([ll]) for a not to be disturbed from its constant value on the inflationary valley of the 
potential. The total energy density of the fields remains smaller than the radiation energy 
density by a factor ej ~ 1.79 x 10 -2 . Moreover, the largest fraction of this energy density is 
stored in oscillations of the ^1,2 fields, whose amplitude is larger by a factor ~ 10 2 than that of 
the <j>x t 2 fields. It is easy then to convince oneself that the constraint of eq. ( pi] ) is comfortably 
satisfied at the onset of inflation. Moreover, as we discuss in the following section, the low- 
frequency oscillations of the <px,2 fields do not generate an effective mass term for a. Only the 
high-frequency ones contribute to this effective mass, and they are much more suppressed at a 
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time tmpi ~ 10 6 . We conclude that inflation sets in as soon as the energy density falls below 
/i 4 , without the requirement of a fine-tuned configuration at the end of the Planck era. 

3 Discussion 

It is possible to obtain an approximate analytical solution for the field evolution during the later 
stages of particle production by the oscillating </> 1)2 fields. At the end of the Planck era the 
energy density is ~ m pl and the "friction" term in the equations of motion is dominated by the 
contribution proportional to 3H. However, after a few oscillations of ^i j2 the Hubble parameter 
becomes smaller than T. The amplitude of the oscillations falls exponentially with time as the a 
field varies very slowly. After a time tmpi ~ 10 2 the amplitude of <^i j2 has dropped sufficiently 
for the evolution equations to be approximated as 

0l,2 + r^ li2 + M|0 1)2 = /i i2 (i) ; (27) 

with M? 2 ~ M? = 2X 2 a 2 . The source terms /i )2 are determined by the last term in the third 

1,2 „/,2> 



line of eq. (|15|): f\ = —KXa^f — ^ 2 )) h = — 2kA<tV ; iV ; 2- The a field is almost constant, while 
the ifti 2 fields oscillate with a characteristic period ~ ir/m ~ 300 m p j, much larger than the 
period ~ tt /mpi of the </>i )2 oscillations. 

An approximate solution of eqs. ([??]) is given by 



M*) = Mt) + Hi,2{t) = '-^r + foe- 2* cos (M^t + a) . (28) 

The term 5<j)i t 2 of this solution describes the decaying high-frequency oscillations, while the 
term 0i 2 determines the variation of 0i i2 induced by the slower evolution of "01,2- It can easily 
be checked that, for 4>X,2 = 4>i,2, the potential becomes independent of a. This implies that 
the effective slope of a during the later stages of the evolution is generated only by the high- 
frequency oscillations 5(f>i^{t). It is their amplitudes that must satisfy the bound of eq. (|lT|). 
Due to the exponential decay of these oscillations until the onset of inflation, their amplitudes 
are smaller than the upper limit allowed by eq. ( |IT| ) by several orders of magnitude. 

A few more remarks are in order with respect to the <r-field evolution before the onset of 
inflation. Its behaviour is not significantly affected by the presence of the gas of ^ particles. 
At early times its evolution is governed by the interaction with the oscillating 0i 2 fields. The 
effective interaction with the thermal gas of scalar ^ particles can be described by the last term 
in the third line of eq. (|i"5|) if the products involving tpi^ fields are replaced by their thermal 
averages. The resulting contribution is zero. The interactions with the gas of fermionic ^ 
particles involve massive $ particles, which are not expected to be abundant. As a result, the 
a field is insensitive to the presence of the fermionic particles. 

The supergravity corrections are not expected to modify the above picture significantly. The 
expansion parameters for these corrections are a 2 /2m Pl , tftl^/'^'mpi, V ; i, 2 /2"7.p;- For our model 
high powers of a 2 /2m Pl are significant, but these are multiplied by powers of /U 2 /m pl near the 
flat direction <j)\ 2 = ipi^ = 0. The most important supergravity correction for our discussion 
arises when the term k 2 (ip 2 + V'l) 2 /4 i n the potential is multiplied by a correction a 2 / (2m' pl ) 
originating in the Kahler potential. A contribution to the effective mass of a is obtained if the 
thermal average of the (ipf + ip 2 ) term is considered. One obtains [21] 



m 2 a (T) = f-?*-, (29) 
m Pl 



7 



with 

The evolution of a as a function of the energy density is determined by the Euler equation plf| 
2 d 2 a t {3 1 ^ dcr / 



0h TT + U - TT — 1— + on I N 2 a = °> ( 31 ) 
\2 1 + tu/ apR 3(l + w) z 



2/ 



with the approximate solution 

a = Gi\ 

.PRi 

Here pm/mpi ~ 10~ 5 , cij ~ 1 correspond to the time when the term of eq. ( |29"|) starts to 
dominate over the mass term ~ A 2 (<^>f + c^ 2 ,) from the potential of eq. (|l5|). The decrease of 
a generated by this supergravity-induced term until the onset of inflation is expected to be 
Aa/mpi ~ 0.2 and should not affect the qualitative features of our study. 

We conclude that the decay of the fields orthogonal to the inflaton into \P particles leads 
to the onset of hybrid inflation without the need of an extremely homogeneous initial field 
configuration. The gas of particles produced before the inflationary era is diluted by the 
subsequent expansion and their number density is negligible today. 

It should be pointed out, however, that a certain degree of homogeneity is still necessary. 
The field evolution must start closer to the flat direction of the potential than to its minima. 
According to the discussion at the end of the first section, this requirement must be satisfied 
over distances far beyond the initial Hubble length. This means that some initial homogeneity 
must exist and extend far beyond the typical size of regions that can be considered causally 
connected. This problem can be resolved if inflation starts at an energy scale near mpi, as for 
example in the models of two-stage inflation of refs. [20, pjf[ . 
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Fig. 1: The evolution of a, 4>i, for a theory described by the potential of eq. (|j) with A = 0.3, 
n/mpi = 10 -3 . The initial values of the fields have been chosen a/mpi = 3, cpi/mpi = 0.6, 
fa/mpi = 0.5. Dimensionful quantities are given in units of mpi. 
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Fig. 2: The evolution of a, (pi, for a theory described by the potential of eq. ( |l5| ) with 
k = 0.5, A = 0.3, n/mpi = 10 -3 , m/mpi = 10~ 2 . The initial values of the fields have been 
chosen a/mpi = 3, (f>i/mpi = 0.6, fo/mpi = 0.5, ipi/mpi = 0.1, ipi/mpi = 0.2. Dimensionful 
quantities are given in units of mpi . 



t 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




